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^\ ABSTRACT 
O 

\-{ ' We consider the spin angular momentum evolution of the accreting compo- 

I nents of Algol-type binary stars. In wider Algols the accretion is through a 

disc so that the accreted material can transfer enough angular momentum to 

^ . the gainer that material at its equator should be spinning at break-up. We 

(N 



demonstrate that even a small amount of mass transfer, much less than re- 



' quired to produce today's mass ratios, transfers enough angular momentum 

to spin the gainer up to this critical rotation velocity. However the accretors 



in these systems have spins typically between 10 and 40 per cent of the crit- 
ical rate. So some mechanism for angular momentum loss from the gainers 
is required. Unlike solar type chromospherically active stars, with enhanced 
magnetic activity which leads to angular momentum and mass loss, the gain- 
ers in classical Algols have radiative envelopes. We further find that normal 
radiative tides are far too weak to account for the necessary angular momen- 
tum loss. Thus enhanced mass loss in a stellar wind seems to be required to 
spin down the gainers in classical Algol systems. We consider generation of 
magnetic fields in the radiative atmospheres in a differentially rotating star 
and the possibility of angular momentum loss driven by strong stellar winds 
in the intermediate mass stars, such as the primaries of the Algols. Differen- 
tial rotation, induced by the accretion itself, may produce such winds which 
carry away enough angular momentum to reduce their rotational velocities to 
the today's observed values. We apply this model to two systems with initial 
periods of 5d, one with initial masses 5 and 3Mq and the other with 3.2 
and 2M0. Our calculations show that, if the mass outflow rate in the stellar 
wind is about 10 per cent of the accretion rate and the dipole magnetic field 
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is stronger than about 1 kG, the spin rate of the gainer is reduced to below 
break-up velocity even in the fast phase of mass transfer. Larger mass loss 
is needed for smaller magnetic fields. The slow rotation of the gainers in the 
classical Algol systems is explained by a balance between the spin-up by mass 
accretion and spin-down by a stellar wind linked to a magnetic field. 

Key words: binaries: close, stars: evolution, stars: magnetic fields 



1 INTRODUCTION 



Pols et al 



19951). There 



Evolution of single stars is now well modelled (see for example 
remain concerns with mass loss, rotation and convection but appropriate and successful 
empirical treatments exist. Evolution of a binary star has several additional complications 
associate d with interaction between the components. Since solving the mystery of Algol 



systems ( IHoyle 



1955 



Crawford 



19551 ). the prototype of semi-detached Algol- type binary stars 



with one evolved and one main-sequence component, we realize that there are some stages 
of evolution when interaction between the components is unavoidable. We must therefore 
take into account, in our calculations, the mass transfer and mass loss together with any 
angular momentum and magnetic interaction between the components, at least in some 
critical phases, to fully understand evolution of a binary system. 

Over the last few decades, the evolution of Algols has been modelled with well defined 
approximations such as conservation of total mass and total orbital angular momentum. 
The effect of mass transfer on the structure of both stars can be modelled reasonably well. 
The angular momentum transfer during mass exchange, however, is not well understood. As 
we shall see in section 1271] there are some episodes of mass transfer in Algols when accretion 
discs or disc-like structures form around the mass gainer. 

Current approximations of binary star evolution do not adequately explain the spin angu- 
lar momentum of the mass-gaining components because the high specific angular momentum 
of the disc material should easily spin these stars up to their critical break-up rotational 
velocities in less than the time needed to reverse the mass ratio of system and enter the 
Algol phase. Here we discuss formation of discs in classical Algol systems and consider the 
spin angular momentum evolution of mass accreting components, taking into account discs, 
tides and magnetic stellar winds. We demonstrate that tidal effects play a minor role in the 
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Table 1. The absolute parameters of Algol primary components dlbanoglu et al.ll2006l V For each star, columns 2 — 8 are the 
orbital period, mass ratio, masses, radii and inclination. Then Vnyn would be the equatorial velocity of the mass gainer (star 1) 
if it were synchronous while Dcq sini is the measured projected velocity and F = Doq/fsyn. 



Name 


P/d 




9 


Mi/Mq 


M2/MQ 


Ri/Rq 


-R2/R0 


i 


fsyn sini 


Dcq sin i 


F 


Ref. 


TW And 


4.12 





210 


1.68 


0.32 


2.19 


3.37 


87 


27 


32 


1.19 


5 


KO Aql 


2.86 





217 


2.53 


0.55 


1.74 


3.34 


78 


30 


41 


1.37 


4 


IM Aur 


1.25 





311 


2.24 


0.76 


2.57 


1.74 


75 


101 


135 


1.34 


1 


R CMa 


1.14 





170 


1.07 


0.17 


1.50 


1.15 


80 


76 


98 


1.29 


1 


S Cnc 


9.48 





090 


2.51 


0.23 


2.15 


5.25 


83 


12 


174 


14.5 


5 


RZ Cas 


1.20 





351 


2.10 


0.74 


1.67 


1.94 


83 


62 


87 


1.40 


3 


TV Cas 


1.81 





470 


3.78 


1.53 


3.15 


3.29 


79 


90 


79 


0.88 


3 


U Cep 


2.49 





550 


3.57 


1.86 


2.41 


4.40 


88 


56 


437 


7.80 


3 


RS Cep 


12.4 





145 


2.83 


0.41 


2.65 


7.63 


87 


11 


170 


15.4 


1 


XX Cep 


2.34 





150 


2.03 


0.33 


2.12 


2.25 


85 


46 


47 


1.02 


2 


U CrB 


3.45 





289 


4.74 


1.46 


2.79 


4.83 


82 


48 


60 


1.25 


1 


SW Cyg 


4.57 





190 


2.50 


0.50 


2.60 


4.30 


83 


22 


196 


8.91 


3 


WW Cyg 


3.32 





310 


2.10 


0.60 


2.00 


7.00 


89 


31 


41 


1.32 


2 


TW Dra 


2.81 





470 


1.70 


0.80 


2.40 


3.40 


86 


43 


37 


0.86 


2 


AI Dra 


1.20 





429 


2.86 


1.34 


2.17 


2.42 


78 


83 


85 


1.02 


1 


S Equ 


3.44 





130 


3.24 


0.42 


2.74 


3.24 


87 


40 


52 


1.30 


4 


AS Eri 


2.66 





110 


1.92 


0.21 


1.57 


2.19 


80 


29 


36 


1.02 


5 


RX Gem 


12.2 





254 


4.40 


0.80 


4.80 


7.00 


85 


20 


157 


7.85 


2 


RY Gem 


9.30 





193 


2.04 


0.39 


2.38 


6.19 


83 


13 


70 


5.38 


5 


AD Her 


9.77 





350 


2.90 


0.90 


2.60 


7.70 


84 


13 


143 


11.0 


2 


TT Hya 


6.95 





224 


2.63 


0.59 


1.95 


5.87 


84 


15 


164 


10.9 


2 


5 Lib 


2.33 





345 


4.70 


1.70 


4.12 


3.88 


81 


89 


68 


0.76 


1 


AU Men 


11.1 





199 


5.93 


1.18 


5.28 


10.04 


79 


24 


124 


5.17 


5 


TU Mon 


5.09 





210 


12.60 


2.70 


5.60 


7.10 


89 


56 


153 


2.73 


2 


AT Peg 


1.15 





484 


2.50 


1.21 


1.91 


2.11 


76 


80 


82 


1.02 


1 


/3 Per 


2.87 





217 


3.70 


0.81 


2.74 


3.60 


82 


51 


52 


1.02 


3 


RW Per 


14.2 





150 


2.56 


0.38 


2.80 


7.30 


81 


10 


161 


16.1 


3 


RY Per 


6.86 





271 


6.24 


1.69 


4.06 


8.10 


83 


30 


213 


7.10 


1 


Y Psc 


3.77 





250 


2.80 


0.70 


3.06 


3.98 


87 


37 


38 


1.03 


3 


RZ Set 


15.2 





216 


5.50 


1.50 


11.00 


14.00 


83 


36 


222 


6.17 


3 


V356 Sgr 


8.89 





380 


12.20 


4.70 


8.50 


15.40 


85 


48 


212 


4.42 


1 


V505 Sgr 


1.18 





520 


2.68 


1.23 


2.24 


2.17 


80 


85 


101 


1.19 


1 


U Sge 


3.38 





370 


4.45 


1.65 


4.00 


5.00 


90 


60 


76 


1.27 


3 


A Tau 


3.95 





263 


7.18 


1.89 


6.40 


5.30 


76 


80 


88 


1.10 


1 


TX Uma 


3.06 





248 


4.76 


1.18 


2.83 


4.24 


82 


46 


63 


1.37 


2 


Z Vul 


2.45 





430 


5.40 


2.30 


4.30 


4.50 


89 


89 


135 


1.52 


1 



References: (1) Ivan Hamme fc WilsonI lll990h (2) lEtzel fc OlsonI l|l993l ').f3'l iMukheriee. Peters, fc WilsonI l|l996l '). (4) 
ISovdugan et al ] l|2003) and (5) iGlazunova et al.l (120081 ). 

removal of excess angular momentum from the gainers and rely on a magnetically locked 
stellar wind to do this. 



2 OBSERVATIONS AND MOTIVATION 

Recently 



Ibanoglu et al. 



( 120061 ) compiled and analysed 61 Algols for which the fundamental 
parameters are well known. In Fig. [1] we show locations of the primary and secondary 
components in the Hertzsprung- Russell diagram (HRD). We use the observers' convention 
and refer to the brighter, hotter and currently more massive star as the primary component 
with mass Mi and the redder, mass losing component as the secondary with mass M2. 
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Figure 1. Locations of the components of the well known Algol systems in a Hertzsprung-Russell diagram. Dots and circles 
arc the primary and secondary components, respectively. The zero-age main sequence ( ZAMS, con tinuous), terminal-age main 
sequence (TAMS, dashed) and base of the giant branch (BGB, dotted) are as found bv lPols et ahl |[T99aD . 



Ibanoglu et al.l (120061 ) arrived at several interesting observations concerning the relation 
between orbital angular momentum and mass. 

(i) Semi-detached binaries (SDBs) with mass ratios q = M2/M1 > 0.3 and orbital periods 
of P > 5d have almost the same angular momentum as detached binaries (DBs). However 
the SDBs with short periods have lower angular momentum even when they have the same 
mass ratios. 

(ii) The orbital angular momenta of SDBs with periods P < 5 d and P > 5 d are 45 and 
25 per cent smaller, respectively, than those of DBs of a total mass of about SMq. 

(iii) The secondaries of SDBs with orbital periods longer than 5 d have angular momenta 
twice that of secondaries with the same mass but with a period shorter than 5 d. 

(iv) The specific angular momenta of systems with P > 5 d are about 24 per cent larger 
than those of the systems with P < 5 d for primary components of the same mass. More 
extremely the specific angular momenta of the longer period systems are 65 per cent greater 
than those of the shorter period systems with the same mass secondary star. 

These results suggest that different mechanisms govern angular momentum evolution 
for short and long period SDBs. Absolute parameters and projected rotational velocities 
measured with great accuracy are given for some Algols in Table [1] The spins of the primary 
stars in SDBs show a marked distinction between short and long period systems at an orbital 
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Figure 2. (a) The ratios of spin angular to Keplcrian angular velocities for the mass accreting primaries of the semi-detached 
binaries are plotted against their orbital periods. The solid line is the synchronous value at that period, (b) The ratios of the 



observed equatorial to the computed synchronous rotational velocities, F : 



q/vsyn, for the Same stars are plotted against the 



mass of the gainer. Open circles and solid dots represent the SDBs with periods P < 5 d and P > 5 d, respectively. There is a 
marked difference above and below P = 5d. At shorter periods the accretors are close to synchronous, while at longer periods 
they tend to spin somewhat faster. 



period of about 5d ( jEggletonll2006l ) as we show in Fig. [21 This indicates that the observed 
rotational velocities of the primaries are far from synchronized with the orbit in longer 
period systems, P > 5 d. The main discriminator of the high rotational velocity of a gainer 
in a classical Algol appears to be its orbital period, independent of its mass. However, if we 
consider only long period Algols, the spin rate appears to decrease as the mass of the gainer 
increases. It should be noted that the gainer in U Cep seems to rotate faster than other short- 
period Algols despite its 2.5 d orbital period. Photometric and spectroscopic observations 
indicate that this system shows transie nt disc because the eclipse durations vary from time 



to time (IGimenez 



1996 



Manzoori 



20081 ). It shows variations in both orbital period and total 



luminosity consistent with mass transfer and convective activity. This slightly anomalous 
behaviour is probably due to the mass ratio in U Cep being very clos e to the critical mass 
ratio at which mass transfer proceeds dynamically (ITout fc Hall! 1 19911 ). 

The above results led us to reconsider tidal interaction and angular momentum transfer 



in systems in which mass transfer is still occurring. A number o 



mechanisms (jPacket 



1981 



Eggleton 



2000 



angular momentum loss 



Chen. Li &: Qianll2006l ) have been suggested to 



explain the evolution of Algols but none is entirely satisfactory. In any case, it is well estab- 
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lished that accretion discs can be formed in phases of evolution when the relative radius of 
the mass accreting star is small enough. 



2.1 Accretion Discs 



Classical Algols are semi-detached interacting eclipsing binary stars in which the less massive, 
evolved secondary component (spectral type F or later G and luminosity class of giant or 
sub-giant) has expanded enough to fill its Roche lobe. These less massive cool secondaries 
are transferring material through a gas stream on to a B or A spectral-type main-sequence 
primary component. In the long period, P > 5d, Algols the mass gaining components are 
small enough, relative to the binary separation, that mass transfer takes place through an 
accretion disc. The in-falling mater ial has too much aii gular momentum for the stream 



to directly impact on the accretor. 



Lubow fc Shul (119751 ) examined the condition for the 



formation of discs in semi-detached systems by modelling the str eam as a ballistic flo w from 
the inner Lagrangian point. As seen in Fig. |3^ (similar to fig. 4 of lLubow fc Shulll975l ). if the 



minimum distance of the stream from centre of the gainer aromm is smaller than the radius of 
detached component {Ri) the transferring material can impact directly on its surface. This 
impact leads to the formation of variable accretion structures. Otherwise, when awram > Ri, 
the mass flow misses the star and collides with itself at a larger radius. Radial motion is 
dissipated and the resulting ring of material spreads viscously to form a permanent accretion 
disc of radius accd- If the accretor has a radius between azu^i^ and aw^ a disc could exist 
because it could intercept the stream before it impacts the star. In such a case the disc 
could be transient. If Ri > aw^ then the stream must impact the star directly. Expressing 
these distances relative to the separation a, we can predict the presence of discs in semi- 
detached systems. Both w^i^ and Wd are functions only of mass ratio q. The majority of 
Algols with P > (4 — 5) d are within the region (Fig. [3b) where we expect stars to have 
either a permanent or transient disc. 

Indeed the prese nce of a disc around the gain er is easily confirmed with optical spectra 



of the primary star. 



Richards fc Albrightl (119991 ) observed and analysed optical spectra of 



Algols and demonstrated that those with periods P > (4 — 5) d show double-peaked Ha 
emission, characteristic of accretion discs, while those of shorter peri ods show only single 



peake d emission, characteristic of a mass stream structure (fig. 3 of [Richards fc Albright 
199% . 
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x/a q=M2/M, 

Figure 3. (a) The orbit of the mass flow in a semi-detached system with masses 5Mq and 3Mq and an orbital period of 5d. 
The coordinates x and y are scaled to the orbital separation a and are in the the equatorial plane of the binar y, perpend i cular to 
the axis of rotation and in a frame corotating with the system. The mass stream trajectory is calculated with lFIannervl l ll975l Vs 
approximations, (b) The positions of all well known Algol primary stars in a fractional radius against mass-ratio diagram. The 
radii below which a disc must form and below which a disc may form ro^j are indicated. The solid dots are the gainers 

with permanent accretion discs among the long period Algols. 



It is thus well established, both observationally and theoretically, that long period Algols 
develop a permanent or transient accretion disc or disc-like structure during mass transfer. 
In the next section we examine angular momentum transport and loss mechanisms to explain 
the asynchronous rotational velocities of detached components in the presence of discs or 
disc-like structures. 



3 MODELS 

In acc retion disc theory, the in-falling matter first forms an accretion ring. 



Stothers fc Lucy 



(I1972I ) proposed a model in which the accreting star is driven into differential rotation by 
the presence of such a ring or disc around it. Viscous forces transfer most of the angular 
momentum (AM) to the outer edge of the ring while mass falls inwards. The ring spreads 
out. Eventually a disc forms and this allows the matter at the inner edge to fall on to the 
surface of the star. For such a Keplerian disc, the angular velocity of material at radius 



R is given by 



nl^^. (1) 
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where G is Newton's gravitational constant and the M is the mass of the accreting star. 
Hence, the specific angular momentum of accreted material at the surface of the star, of 
radius R, is 



VGMR. 



(2) 



This specific angular momentum of accreted material would be equal to that at the equator 
of the gainer if it were critically rotating at brake up. It is much larger than found in normal 
stars. When accreting at a rate Mace the rate of angular momentum transferred from the 
disc to the star is 



dt 



(3) 



Assuming a negligible change in stellar radius, we can determine the amount of mass 
AM that must be transferred through the disc to spin the star up to its critical angular 
velocity Qy^ from an initial Qo when it had mass Mq. Let the radius of gyration of the star 
be kR so that its total angular momentum is k^MR^Q when spinning rigidly at Q then 

k^ ( . Vln 



AM 



A more precise formula was derived by 



fc2 



Packet 



Jo 



Mn. 



(4) 



(1198 ll ) who took account of the change in 
the mass of the star but this is unnecessary for our purposes because AM is always small. 
For main-sequence stars k"^ ^ 0.1 and varies little. Thus when 0.1 < Qq/Q]^ < 0.4 we find 
0.1 > AM/Mq > 0.06. This is very small when w'e consider that all c lassical Algols have a 
mass ratio of g < 0.7 (mostly q ~ 0.2 according to Ibanoglu et al.l2006 ) which indicates that 
the losers in the classical Algols have transferred more on the order of 1 M© . The shaded area 
in Fig. m shows the amount of material that must be accreted from a disc to spin the star 
up to its critical rotational velocity. Despite having high spin velocities, 0.1 < Q/Q]^ < 0.4, 
observations show that the detached components in most of the Algols do not actually attain 
their critical rotational velocity (Table [H Fig. [2|i). The only exception, wi th the high ratio o f 
Q/Q]^ = 0.72, is RZ Set. This star's radi al velocity curve i s distorted (e.glMcNamaralll957l) . 



It shows emission in outside eclipse ( iMcNamara 



1957 



Hansen &: McNamara 



Olson fc Bell 



19591) and 



19891). These 



its light curve displays distortions due to an accretion stream (|( 
observed phenomena may be taken as the signature of a rapid mass transfer phase. 

In all mechanism is needed to dissipate this excess angular momentum, along 

with associated energy. Here we examine various mechanisms for angular momentum loss 
and compare with the observations. 
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Figure 4. The amount of accreted mass needed to spin a star with initial mass of Mq/Mq up to its critical rate when 
initially f!o = 0.1 (solid line) and 0.4^1; (dashed line). 



3.1 Tidal Forces and Energy Dissipation Mechanisms 



Tidal forces have been quantified since the work of iDarwid (1l879l ) and are, perhaps, the 
most well understood mechanism for angular momentum transfer within binary stars. Tidal 
interactions exchange angular momentum between the orbit and stellar s pins by torq ues. 
They act to synchronize stellar spins with the orbital period. According to IZahnI (l2005l ) the 
time-scale for synchronization tgync is given by 



1 



1 



ts 



-O^^ - J^r-^ - , (5) 

■-sync \L — OJ Idiss -< V ^ / 

where u is orbital angular velocity, a is separation of system, I = k'^MR^ is the moment 
of inertia of the primary star and tdiss is the time-scale for the most effective dissipation 
mechanism. It depends on the stellar type. For stars with convective envelopes the kinetic 
energy of the equilibrium tides is dissipated by turbulent convective eddies. For those with 
radiative envelopes the shortest time-scale is through gravity wave dissipation. It is well 
known that the ene rgy dissipat ion in convective envelopes is much more effective than in 
radiative envelopes (jZahrul2005l ). 

As we shall describe in the details of the binary star evolution models in section 13.21 
we have examined the efficiency of the tidal torque at opposing the spinning up effect of 
accreted material from a disc. Gainers in Algol systems are early type stars, so we take the 



approximation of the dynamical tide with radiative damping given by iHurley. Tout fc Pols 
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(120021 ) to evaluate tdiss for equation [51 In the classical evolution of an Algol the initially more 
massive star evolves more quickly, overfills its Roche lobe and begins a phase of rapid mass 
transfer soon after it departs from the main sequence as it rapidly becomes a giant star. The 
primary star is at this time relatively compact and substantially smaller than its companion 
so that an accretion disc is likely to form. When the accretion is conservative and uniform 
over time, so that the rate of accretion from the disc to the gainer (Mace) is the same as the 
mass- loss rate from the donor (M2), we can calculate the spin angular momentum gained 
during mass transfer. Initially the angular momentum of the gainer is 

Jso = k'^MR^Vlo . (6) 

The net change depends on the competition between the torque of accreted material (equa- 
tion [3]) and the opposing spin-down torque of tides obtained from equation |5l by 

(T) 

where k remains approximately constant. In Fig.O we show the variation in the angular 
velocity of a gainer undergoing mass accretion through a disc with its increasing mass. We 
used a detailed evolution model for a system of stars initially of masses 5 and 3 Mq and 
an orbital period of P = 5 d. We allow the star to rotate faster than its break-up rate for 
illustration only. We find tides are almost incapable of synchronizing the star with the orbit 
because the gainer reaches brake-up velocity after accreting only a small amount of matter 
as indicated by equation HI We found that, to affect the star's spin, the tides would need 
to be stronger by more than a factor of 10''. There is no known physical basis for this. We 
might also hope that a convective core's ability to dissipate energy by tidal forces may have 
an effect but the strong dependence of the synchronization time on {R/a)^ in equation [5] 
means that the contribution of the convective core is too small and may be neglected. Thus 
we conclude that tides are insufficient to synchronize the stellar spin of the gainer with the 
orbit when accretion is through a disc. 

3.2 Magnetic Winds 

Our analysis has shown that tides are insufficient to dissipate the excess angular momentum 
accreted by Algol primaries with discs. So we consider alternatives that are applicable to 
this case. The total angular momentum lost from a star in a wind coupled to a magnetic 
field is equivalent to a ngular momentum ca r ried aAvay by the wind material co-rotating up 



to the Alfven surface (IWeber fc Davis 



1967 



Mestel 



19681 ). So the rate of change of angular 
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Figure 5. Spin evolution of a gainer in a binary with initial masses of 5 + 3Mq and a period of P = 5d in the case when 
total mass and angular momentum are conserved. Calculations were made for standard tides (equation [5] solid line) and tides 
increased by factors of 10^ (dot-dashed line) and 10* (dashed line). The fully synchronous case (dotted line) is plotted for 
comparison. 



momentum of the star owing to the wind is 



(8) 



where Q is the angular velocity of the star and Ra is the Alfven radius at which outflow 
speed equals the local magnetic Alfven speed 

where pA is the density of wind material, Ba is the magnetic field strength at the Alfven 
surface and M„ < is the mass-loss rate. For a spherical outflow 

-Mw 



PA 



AttRIva ■ 



(10) 



The angular momentum loss rate depends on the field structure and flow velocity which 
are not easily determined a priori for the wind. We need to make assumptions about both. 
To model the magnetic field structure in a simple manner, we assume that field strength 
follows a single power law of the form 

Ba = B,iR/RAr, (11) 



wher e n describes the geometry of the stel 



field ( IWeber fc Davis 



1967 



Mestel fc SpruitI 



ar fi eld and n = 3 corresponds to a dipole 



19871 ) while Bs and Ba are the magnetic flux 



densities at the stellar surface and at the Alfven radius, respectively. 

It is usually assumed that the thermal wind velocity is of the order of the escape velocity 
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flTout k Prinde 



1992 



Stepienlll995l ) so that 



2GM 



Rfs 



Combining equations [121 El [10] and [9] we find equation [H] becomes 



(12) 



(13) 



The second and considerable torque arises between st ar and disc by niagneti c interaction. 
The framework of such an interaction was constructed by lGhosh fc Lamb! (119781 ). Some of the 
stellar magnetic dipole flux connects to the accretion disc and transports a ngula r momentum 



betwee n star and disc . We us e the expression given by 
as did 



Armitage fc Clarke! (Il996l ) and assume, 



radius. 



Stepieril (j200d . |2002|), that the radius of magnetosphere is equal to the co-rotation 

(14) 



Rc 



As they pointed out the disc torque does not depend on the mass of the disc so 



'disc 

dt 



3GM' 



(15) 



Stepieh 



where n = BgR^ is the magnetic moment of stellar magnetic field. Contrary to 
( I2OO2I ) who assumes magnetic flux is constant, we assume that the magnetic field strength 
Bs remains constant because the mass of the accreting star increases at a substantial rate 
(Mace ~ lO~^M0yr~^). We seek a solution for various fiducial values of Bs- 

The third torque, the accretion torque discussed in section [HI has the effect of spinning- 
up the mass-accreting star. Now we replace the radius of the star R in equation [3] with 
RcoT because the stellar magnetic field disrupts the disc at i?cor as constant at the value of 
co-rotation radius. 

The wind mass-loss rate remains an unknown parameter. Although there is no a priori 
rate we can set an upper limit. All observed Algols show reversed mass ratio so much of the 
material lost by the donor must be accreted by the gainer. We may write 

Mace ~ /3M2 (16) 

and 

(1-/3) 



M„ < (1 - /3)M2 



(17) 



Matt fc Pudritzl (12005! ) claimed 



where < (3 < 1. For conservative evolution (3 = 1. 
that the mass outflow rate in stellar winds is about 10 per cent of the accretion rate for 
the pre-mainsequence stars, corresponding to /3 ~ 0.9. This is why the classical T Tauri 
stars spin at less than 10 per cent of their breakup velocity. It is the original suggestion by 
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Hartmann fc Staufferl (119891 ) that about one tenth of accreted material lost in a wind can 



remove the accreted angular momentum. 



Such models have a. 



pre-mainse qunce stars flMatt &: Pudritz 



ready been app 



ied to a wide variety accreting stars inclu ding the 



2005 



2008a 



2000 



b|) and the Ap and Be stars (IStepieh 



20021 ). Both Mace and are taken as free parameters in most of these studies. In the 



case of binary evolution Mace depends on the orbital evolution which can be well modelled. 
So, for a given magnetic field, we need only make an estimate of the fraction of mass lost in 
the wind (equation [T7I) to estimate the angular momentum loss. 

Like the tidal torque, which also depends strongly on the separation of binary, the mag- 
netized wind and disc torques are very sensitive to changes in the radius R of the mass 
accreting star so a detailed evolutionary model is required. Here we use binary star mod- 



els ev olved with the Cambridge STARS code whi ch was orig i nally developed by 



Eggleton 



(jl97ll ). The physics was systematically updated by lPols et al.l (119951 ) and the code has been 
modified so that it can evolve both components of a binary system sir nultaneously together 
with the effects of mass and angular momentum loss and transfer (IStancliffe fc Eldridge 



20091). 



We have made various binary evolution models for two different initial systems, each 
with a period of 5 d and with masses 5 + 3 Mq and 3.2 + 2 Mq. Such systems are believed to 
be the progenitors of Algol-type systems. We evolved both with /3 = 0.9, 0.7, 0.5 and 0.1. 
The code does not directly take into account the presence of a disc so we tabulate the output 
of the one dimensional structure of each component, the mass-transfer rate, accretion-rate 
and orbital parameters at every time step. From this we calculate the torque of the magnetic 
wind, the disc and the accretion and hence the resulting rotation rate for a given magnetic 
field strength. Because the variation in radius of the gainer is taken directly from the tables it 
is an approximation. It does not take account of possible rapid rotation. Typically the radius 
of an early-typ e non-rotating star depends on the mass as ~ M"'^'' at solar metallicity 
(IMaederl |2009| ). A rapidly rotating star becomes oblate. Its polar radius becomes smaller 



than that of the equatorial radius so that the equatorial radius, when rotating solidly at 



critical 



velocity is increased by a factor of 1.5 with respect to a non- rotating star of the same 



mass ( lEkstrom et al.l 120081 ). However when Q/Qk ~ 0.8 this factor is only 1.1 and can be 



neglected when compared to the increase in radius as matter accretes. 

To calculate the spin evolution of the mass accreting component we assume a dipolar 
(n = 3) magnetic field and that a magnetized wind is launched as soon as mass transfer 
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Figure 6. Evolution of a binary star with initial mass of 5 and 3Mq and orbital period of 5d with /3 = 0.9 and Bs = 1.5 kG. In 
(a) the dashed line is the track of the initially more massive donor and the solid line that of the gainer in an Hertzsprung-Russell 
diagram. Panels (b), (c) and (d) show the accretion rate, radius and f2/f2fc of the gainer as a function of time since the onset 
of Roche lobe overflow (RLOF). 



begins. We first obtain an evolutionary model of the given binary and each value of /3. Then 
for various field strengths of Bg for each /3 we examine carefully how the angular velocity 
of the star changes as its mass increases. If the star never reaches its critical velocity the 
computation is ended. If, on the other hand, the star reaches its critical rotation rate we 
interrupt the mass accretion but allow the magnetic wind and disc torques to lower the 
angular velocity to less than the critical rate. This leads to an extra mass loss and a time 
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lag in the evolution of the system. When the gainer loses enough angular momentum that 
Q/Qk < 1 mass accretion resumes. 

In Fig. |6]we show the evolution of the system with initial masses 5 and 3 Mq and orbital 
period of 5 d. The 5 Mq star evolves off the main sequence when it exhausts hydrogen in 
its core and moves towards the Hertzsprung gap where it fills its Roche lobe and rapidly 
transfers most of its envelope to its less massive companion. It becomes established as a 
red giant when its mass has fallen to 0.67 M0. Our models show that this final mass is 
almost independent of f3 (0.675 Mq for (3 = 1 falling to 0.665 for /3 = 0.5). Meanwhile, 
the gainer moves up the main sequence as mass transfer proceeds. In Fig. |6^ the dashed 
line is the evolutionary track of the 5 Mq and the solid line that of the 3 Mq star. The 
mass accretion rate Mace for P = 0.9 is plotted against time in Fig. [Hb- It is larger during 
the early stages and slows sharply after 10^ yr as the loser reaches the giant branch. The 
gainer reaches hydrostatic equilibrium on a dynamical time-scale and then thermodynamic 
equilibrium on the Kelvin-Helmholtz time-scale. The variation of the accreting star's radius 
over the same time interval is plotted in Fig. ^ and its angular velocity in Fig. |6]i when the 
magnetic field strength is = 1.5 kG. Recall that the variation of angular velocity is caused 
by the interaction between the torques of the wind, the accretion and the disc together with 
radius changes of the star. After the phase of rapid mass transfer (g ^ 0.2) the slowing 
torques dominate. The location of the donor in the H-R diagram (see Fig. [1]) shows that 
it is still at its minimum luminosity reached following mass transfer. For this system, this 
phase corresponds to an effective temperature of log^Q(Tcff/K) = 3.66 and luminosity of 
logio(L/LQ) = 2.25. 

The angular velocities computed for the same initial system but varying /3 and are 
shown in Fig. [71 As expected a larger Bg creates larger torques for the same /3 and hence 
smaller angular velocities fi/i^fc at a similar point in the evolution. For the same B^ but 
smaller /3 there is more mass lost in the wind and again smaller Q/Qk but then the evolution 
of the system is highly non-conservative. The most important result we find here is that 
when Bg > 1 kG the gainer does not reach its critical spin rate if it loses 10 per cent of the 
mass transferred from its companion. In such classical Algol as observed today is 

produced. For decreasing /3 the final mass ratio is larger for the same initial mass. 

We also applied our models to a system with initial masses of 3.2 + 2Mq and period 



P = 5d. This is typical of the progenitors considered for Algol (/3 Per) itself (iTout et al 



19971 ). The models are shown in Figs. |8]and|9l 
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Figure 7. Variation of f2/f2fe, the fractional angular velocity of a gainer of initial mass 3Mq, with time elapsed since the onset 
of RLOF for (a) /3 = 0.9, (b) 13 = 0.7, (c) 13 = 0.5 and (d) 13 = 0.1 with Bs = 1, 2, 3, 4 and 5kG, from top to bottom in each 
panel. 



Although there is evidence that magnetic winds brake the spin of stars, the underlying 
physics that produces the magnetic fields is not yet fully understood. There is a distinction 
between the magnetic properties of stars on the upper and lower main sequence. In the 
case of the Sun and solar-like stars the combination of a convective envelope and rotation is 
thought to generate the observable magnetic fields by an a—Q dynamo, the strength of which 
depends upon the stellar angular velocity. However, in contrast, only a minority of early-type 



stars, such as the Ap and Bp stars have observable m agnetic fields 
been made to explain these phenomena (for example lAlecian et al. 



Vlany attempts have 



2007 



Ekstrom et al. 
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Figure 8. Similar to Fig. [6] but evolution of a binary star with initial masses of 3.2 and 2M0 and orbital period of 5d with 
/3 = 0.9 and Bs = 1.5 kG. In (a) the dashed line is the track of the initially more massive donor and the solid line that of the 
gainer in an Hertzsprung-Russell diagram. Panels (b), (c) and (d) show the accretion rate, radius and Q/Qk of the gainer as a 
function of time since the onset of RLOF. 



2008 



Ferrario et al. 



20091 ). The primaries in most Algols, especially those with discs, are 



late B or early A type stars and have radiative envelopes. Stellar models confirm that the 
envelope is convectively stable in massive stars. Therefore solar-like concepts of dynamo 
operation are not directly applicable to these stars with radiative envelopes. Furthermore, 
there is no observational evidence of the presence of active regions on the surfaces of A and B 
stars. The origin of any magnetic dynamo which feeds energy into winds in early type stars 
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Figure 9. Similar to in Fig.[7]but for initial masses of 3.2 + 2M0 and period P = 5d, /3 = 0.9, 0.7, 0.5, and 0.1 and Bs = 0.5, 
1, 2, 3 and 4kG. 



is unclear. Similar problems also exist 



within the radiative layers of the Sun (jChaplin et al. 



or the explanation of the solid-body rotation profile 



2001) 



While there is a growing body of observational evidence (jPonati et al 



2006 



MacDonald fc MuUan 



20041 ) that early type stars have measurable magnetic fields, no attempt has led to a com- 
prehensive theory which explains both the origin of and the sustai ning mechanis m of the 
magnetic field. Two leading possibilities are the fossil field theory ( ICowlinglll945l ) and the 
dynamo. Existing theories of dynamo generated fields focus on different assumptions but 
propose that dynamos ope rate either in the fully conve ctive cores or in the differentially 
rotating radiative envelope. iMacDonald fc MuUanl (120041 ) report serious difficulties with the 
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core dynamo hypothesis suggested by ICharbonneau fc MacGregoii (120011 ) because the core 
must create magnetic fields that are much stronger than equipartition values in order to 
reach to the surface. 



Spruit 



(Il999[ ) argues that a principle ingredient of dynamo operation in a star is rapid ro- 
tation. This is certainly available in massive stars. While the driving mechanism of magnetic 
field generation is based on hydrodynamical instabilities in solar type stars, similar insta- 



bilities in the radiative stars may be pr oduced by the magnetic fi e 



are the Parker and Tayler instabilities ( iParker 



197 



Tayler 



197: 



ds themselve s. Examples 



SpruitI ( 120021 ) formulated 



a dynamo mechanism for a differentially rotating star in which Tayler instability replaces 
the role of convection in closing the field amplification loop. The azimuthal component of 
the magnetic field grows until it reaches a critical strength by winding up an initially very 
small radial component in a non-convective zone. At that point the Tayler instability op- 
erates on a very short time-scale, of the order of Alfven crossing time, to regenerate the 
radial field. The result is a predominantly toroidal fie ld an d a closed dynamo loop. Recent 
numerical simulations by iBraithwaite fc SpruitI (l2004j ) and iBraithwaite fc Nordlundl ( 120061 ) 
showed that a stable configuration can be reached within a non-convective star on an Alfven 
time-scale from an arbitrary initial magnetic field. The Spruit- Tayler mecha nism is not the 
only one that operates in non-convective material. iBalbus Sz Hawleyl ( 119941 ) have proposed 
a magneto-rotational instability which works well for accretion discs without convection. 

Fig. [To] shows the internal structure of the gainer during the mass transfer for our system 
with initial masses of 5 and 3 Mq. The gainer has a substantial convective core and an almost 
consta nt mas s radia tive envelope throughout the mass accretion. 

As ISpruitI (120021 ) pointed out, differential rotation in single stars provides a finite amount 



of energy for generating magnetic fields so we might expect that the dynamo to cease once 
the rotation profile has been smoothed out. In the case of our semi-detached binary stars 
with discs, the energy in differential rotation can be continuously supplied from the disc 
material which has high specific angular momentum relative to the star as we discussed in 
section [31 Hence, if a Spruit-Tayler type dynamo operates in this case, magnetic fields are 
regenerated as long as mass transfer continues. According to our assumption that the wind 
is accelerated to the stellar escape velocity (equation [12]) and leaves the system at the Alfven 
radius, the rate of energy required to be fed into the wind is given by 

GMM^ 



Ra 



(18) 
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Figure 10. Radiative and convcctivc zones with mass m in the gainer in our system of initially 5 + 3Mq, orbital period 5d 
and (9 = 0.9. Convective (cross shaded) and non-convective (horizontal shaded) regions of the gainer during the mass transfer 
are shown with respect to the time elapsed since the onset of RLOF. The envelope of the gainer is always convectively stable. 

This can be provided via the dynamo from the highly energetic disc material which con- 
tinuously supplies energy to the star. Roughly speaking, AQ ^ RdQ/dr is the change of 
angular velocity between centre and outer edge of the gainer or across its su rface from pole 



to equator. The shear energy which can be tapped is (iTout fc Pringldll995l ) 

I „ / (t/I//" \ / /\\/. \ / \/. \ ^ 

-B.h » 



1 /GM^ 

2 \ R 



n 



n 



(19) 



If AQ ^ Q we can evaluate Esh/L„ to estimate the time over which the shear energy could 
sustain the wind, even if it were not replenished, to be 10^ — 10^ yr. This is similar to the 
timescale over which mass transfer takes place so the shear energy is a good candidate for 
supplying energy to the dynamo during the mass transfer phase. The quantity rdyn = Esh/L^ 
du ring the mass transfer is plotted against the time since the onset of RLOF in Fig. [TTl 



Tout &: Pringld (1l992l ) developed a model based on the idea that the dynamo process 



creates not only magnetic field also leads a continual expulsion of magnetic flux from the 
star. This flux, they assume, provides the mechanical output of energy that drives stellar 
wind. They developed schematic dynamo equations and give a rough expression for this rate 
as 

4 



^0 -l4 ^3 

471 3 



(20) 



where is the decay time-scale of the toroidal field component B^. It is assume d to be equal 
to Tg, the growth time-scale on which the dynamo reaches equilibrium, given by (jSpruitll2002l . 
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t-tRLOF/y 

Figure 11. Variation of the dynamo timescale r^yn = ^sh/^w with time elapsed since the onset of RLOF. The shear energy 
seems to be sufficient to sustain the required wind for 10^ — 10^ years at any stage. 



Eq 17) 



Combining equations [18] and 



and using equation [21] with equation 19 of 



(21) 



Spruit 



(12002), 



we can estimate the mass-loss rate in the wind for mean values of r ^ R and p ^ p, 

2 / X 1/2 / 7^ \ 1/2 



M 



n'J [nJ [rmJ 



(22) 



In this formula the quantity is the Brunt- Vaisala frequency, g,. = dlnQ/dlnr is dimen- 
sionless gradient of rotation and K is the thermal diffusivity. This equation is derived under 
the assumption that almost all the flux lost is transferred to the wind. Therefore, equation [221 
is a theoretical upper limit to the mass lost in the wind. We find that such a fully efficient 
Spruit-Tayler dynamo can support a wind up to 0.01 M0 yr^^. In this estimate we use mean 
values of the parameters ~ 10~^s~^, K ^ 10® — lO^cm^ s~^ and assume the differential 
rot ation gradient Qr ~ AQ /Q 

i) have also studied accretion powered stellar winds, both 



Matt fc Pudritd fl2005 



2008a 



analytically and with two dimensional magnetic wind solutions. Although they are interested 
in T Tauri type stars, they claim that their models are independent of dynamo mechanism. 
Based on two-dimensional simulation they found that the field strength outside the star falls 
off as if with n ^ 3.24 which is very similar to the simple dipolar fi eld decay with ^ = 3. We 
compare the values of the Alfven radius we calculate with those of [Matt fc Pudrita (l2008a[ ) 
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Figure 12. Comparison o f our calculated Alfven radii (open circles) with those of lMatt fc Pudritj l l2008al . dashed line) and 
lud-Doula fc Owockil | |2002| . solid line) for a range of values of the wind confinement parameter r\. 



and lud-Doula fc Owockil (l2002l ) in Fig. [121 In this comparison we used the dimensionless 
wind magnetic confinement parameter, r/ which characterizes the ratio of magnetic field 
energy density to the kinet i c ener gy density of the wind. We find a good agreement with 



that of 



ud-Doula fc Owockil (I2QQ21)- However the A 



factor of two smaller than those of 



Matt fc Pudritzl (l2008af l. 



f ven rad ii calculated by us seem to be a 



4 CONCLUSIONS 

It has been known for some time that many long-period Algol systems have accretion discs. 
Accreting material from such a disc should increase the spin rate of the more massive com- 
ponent up t o its break-up speed as soon as even a small fraction of the mass has been 



transferred (Ide Mink. Pols fc Glebbeekl l2007l ). All the classical Algols have a less massive 



evolved and a more massive main-sequence component. Therefore a substantial amount of 
mass from the initially more massive star must be either lost or transferred to its companion. 
The angular velocities of the more massive components in many Algols are measured with 
great accuracy. These measurements show that the gainers rotate somewhat more slowly 
than their break-up rates. Thus there must be a mechanism which removes spin angular 
momentum from the rapidly rotating hot star. We have demonstrated that tidal interaction 
between the components and the orbit is too weak to do this fast enough. We show, however. 
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that magnetic braking, driven by a magnetic d ynamo that is maintained by the accretion 
itself, can remove sufficient angular momentum. iMacDonald fc MuUanl ( 12004J ) hypothesized 
that a dynamo operating in a sheared radiative region can generate flux tubes which are 
able to rise the surface of the star. A shear dynamo in the extended radiative envelope of 
a massive star can therefore serve as an efficient supplier of magnetic flux to the surface of 
the star. We propose a similar self-consistent model in which a dynamo operates in shear- 
unstable material throughout the radiative envelopes of the gainers in Algols with orbital 
periods longer than 4 — 5 d. Such stars can lose angular momentum in magnetized stellar 
winds, leading to non- conservative angular momentum evolution. 

Under this non-conservative evolution, an Algol system continues to evolve with rela- 
tively little mass loss at a rate ~ O.lMacc but corotating to a relatively large Alfven 
radius Ra- When it reaches the classical Algol phase the less massive component is observed 
to be more evolved than the massive one. Spruit-Tayler instabilities seem most likely to be 
responsible for the operation of a magnetic dynamo in massive stars with radiative envelopes 
( iMuUan fc MacDonaldl |2005[ ) . However we need to construct two dimensional models with 
an accretion disc to test how such a dynamo really operates in Algols. Such numerical models 
could also supply more realistic Alfven radii based on the field geometry. Using such detailed 
models we should be able to find an equilibrium spin-ratio which can explain the observa- 
tions in Fig |2l We might also take into account the effec t of w ind induced hydrodynamic 



instabilities suggested by 



Lignieres. Catala. fc Mangeneyl ( 119961 ) which may also effect the 



differential rotation parameter. 

Using a very simple model we computed the spin angular momentum evolution of gainers 
with discs in two systems with an initial orbital period of 5 d but different masses. The results 
show that even a small amount of mass, about 10 per cent of the transferred material, lost 
by the gainer with a magnetic field of 1 kG is sufficient to slow down the star to below its 
break-up velocity. As the magnetic field strength increases the rotational velocity of the star 
decreases for the same amount of mass loss. 
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